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Abstract 

We seek a measure of the neutron density of 208 p]-, fj-om analyses of interme- 
diate energy nucleon elastic scattering. The pertinent model for such analyses 
is based on coordinate space nonlocal optical potentials obtained from model 
nuclear ground state densities. Those potentials give predictions of integral 
observables and of angular distributions which show sensitivity to the neutron 
density. When compared with experiment, and correlated with analyses of 
electron scattering data, the results suggest that 208p]^ ^ neutron skin 
thickness ~ 0.17 fm. 
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I. INTRODUCTION 



Interest in the matter distributions of ^°^Pb, and its neutron density profile particularly, 
is quite topical flj. There is a proposal to measure its neutron root-mean-square (rms) radius 
at the Jefferson Laboratory 0] from an analysis of parity-violating electron scattering data. 
In contrast to proton rms radii that are known to within an accuracy ~ 0.02 fm [Q, neutron 
rms radii are less certain. 

Recently, the neutron rms radius in ^''^Pb was assessed in terms of modern Skyrme- 
Hartree-Fock (SHF) models With the Friedman-Pandharipande (FP) neutron equation 
of state [H| as a constraint, the neutron rms radius in ^°^Pb was expected to be 0.16±0.02 fm 



larger than the proton value. Previous estimates of this neutron skin, S = y (r^) — y {Tp 

ranged from 0.1 fm to 0.3 fm 0,^; the lower values favored in general by SHF models, while 
relativistic mean field models predict values closer to 0.3 fm 0. Knowledge of the skin 
thickness in ^°^Pb then is a good constraint upon such models of structure ||^. 

The planned parity-violating electron scattering experiment will only provide infor- 
mation about the neutron rms radius itself. We seek further information and address the 
question of whether analyses of nucleon scattering data establish a measure of the neutron 
density distribution in ^''^Pb. Hadron scattering data have been analyzed previously to de- 
duce the neutron skin thickness in ^°^Pb. The ratio of tt"^ and tt" reaction cross sections gave 
S = O.OibO.l fm |J7[. Analyses of elastic proton scattering data at 0.8 GeV gave 0.14±0.04 fm 
ip. However, a review of the analysis of proton scattering from ^°Ca ||^ gave values in the 
range of —0.4 fm to —0.2 fm for 5* in that nucleus which are systematically smaller than all 



theoretical models which give —0.05 fm |T0[, and suggests that there is a systematic problem 
in the analysis using phenomenological models of high energy proton scattering data which 
would affect the extracted S values at the level of ~ 0.2 fm. A more recent analysis of 
650 MeV proton scattering data gave S = 0.20 ± 0.04 fm for ^osp]^ while that reaction 
model gave a result for ^°Ca which is consistent with theoretical predictions. The excitation 
of the isovector giant dipole resonance in ^osp]^ inelastic alpha scattering ||12[ was used 
to deduce S = 0.19 ±0.09 fm. These previous analyses, based on phenomenological models, 
produce a range of results for ^osp^^ ^^^^ suggest that there may be small but systematic 
errors in those strong interaction models for hadron scattering limiting the accuracy in the 
extraction of S" to ~ 0.2 fm. 

Our approach is based on coordinate space nonlocal optical potentials generated by a 
full folding of realistic effective nucleon-nucleon (NN) interactions with ground state density 
matrices (termed densities hereafter) of ^°^Pb. This allows us to distinguish between various 
model structures, including those of the SHF type proposed by Brown [Q], even if two may 
have the same rms radii. As a calibration of the use of SHF models we consider the elastic 
scattering from ^°Ca as well. Further, as the effective A^A^ interaction is dominated by the 
isoscalar ^Si channel |]13[, proton scattering predominantly will probe the neutron density 
and vice-versa and so we consider both proton and neutron elastic scattering at a given 
energy seeking as complete a map as possible of the nuclear matter distributions. 
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II. NUCLEAR MODELS AND THE MICROSCOPIC OPTICAL POTENTIAL 



The model with which predictions of nucleon-nucleus {NA) scattering observables are 



made has been given in detail in a recent review |jT3[. Use of the complex, nonlocal, NA 
optical potentials defined by that model prescription, without localization of the exchange 
amplitudes, gave predictions of differential cross sections and spin observables that are in 
good agreement with data from many nuclei (^He to ^^^U) and for a wide range of energies 
(40 to 300 MeV). Crucial to that success was the use of effective NN interactions built upon 
A^A^ g matrices; solutions of Brueckner-Bethe-Goldstone equations for realistic starting (free) 
A^A^ interactions. The NA optical potentials result from folding those effective interactions 
with the densities of the target nucleus. That folding includes the antisymmetrization of the 
projectile-nucleus wave functions and therefore exchange (knock-out) amplitudes are treated 
explicitly. Consequently the NA potentials inherently are non-local. The optical potentials 
that result from that process we term for brevity as (^-folding potentials. 

Recently, this approach was applied successfully to make predictions of the integral ob- 
servables of nucleon elastic scattering [T^. Thus, the use of the ^f-folding optical potentials 
give good predictions to both angular-dependent and integral observables; a result not guar- 
anteed with the more common phenomeno logical approaches. Of import, however, is that 
the level of agreement with data in the (yf-folding model depends on the quality of the un- 
derlying model of structure. We seek to use that dependence as a sensitive evaluation of the 
densities considered. 

Our theoretical density distributions are based upon the Skyrme Hartree-Fock model 
for 208p|^ with a spherical closed-shell configuration. It was shown in that there is 
a combination of parameters in the Skyrme Hamiltonian (dominated by the X3 parameter) 
which has a strong influence on the neutron skin thickness but which are not well determined 
by the data on binding energies and charge radii. This combination can be correlated with 
the pressure in the neutron equation of state at normal nuclear density It can also be 
related to the surface symmetry energy |T^. The Skyrme interaction SKX |]T5| was obtained 
with a constraint on the neutron equation of state provided by the FP model which 
constrains the neutron skin in ^°^Pb to be S" = 0.16 ± 0.02 fm. The results obtained with 
the SKX model will be denoted herein by SHFl. If one does not allow for any constraint 
from a model for the neutron equation of state then a much larger range of S is allowed. In 
particular, we use a model which is constrained by the same nuclear properties used for SKX 
but with a value of the parameters which gives 5* = 0.25 fm for ^°^Pb. This we denote as 
SHF2. The densities for ^osp]^ obtained from these models, as well as the others considered 
herein, are shown in Fig. |I]. 

An initial test of these interactions is provided by the elastic electron scattering data 
which yields information on the charge density which, in turn, gives information on the 
proton density. SKX appears to give an excellent reproduction of the charge-density dis- 
tribution (Fig. 11 of [|l3). However, there is some model dependence in the extraction of 
the charge density from electron scattering data. It is better to compare to a representation 
which is more closely associated with the actual data - this is the plane-wave transform of 
the charge density shown in Fig. The experimental form factor is obtained from the charge 
density distribution given in ||^. Experiment is compared with the SHFl and SHF2 model 
results showing a disagreement with data which systematically increases with momentum 
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transfer. The main feature of the distribution which affects the high-g behavior is the sur- 



face thickness. As discussed in ||15[, the SKX interaction appears to have a surface thickness 
which is a httle sharper than that determined from experiment. We have thus looked at 
other Skyrme interactions in terms of the data in Fig. |] and find the older SKM* interaction 
||T6| is much better than others in this regard. SKM* appears to achieve this improvement 
by a decrease in the power of the density, p", associated with the density-dependent part 
of the interaction from its value a = 1/2 for SKX to a = 1/6 for SKM*. Coinci dentally, 
SKM* predicts a neutron skin of 0.17 fm which is essentially the same as that predicted by 
SKX (SHFl model). This is obtained mainly because the parameter was set to zero by 
default, since it is not well determined by nuclear data. 

The value of a is also associated with the incompressibility coefficient K for infinite 
nuclear matter which ranges from K = 270 MeV for SKX to = 217 MeV for SKM*. It 



was found in that a = 1/2 gave the best overall fit with the data set considered and 
that when a is decreased (as in the SKXm interaction, which has a = 1/3) the overall 
increased mainly because of an increase in the contribution from the single-particle energies. 



It was also argued in |15] that an improvement in the surface properties may require an 



additional parameter in the Skyrme Hamiltonian associated with the next order d-wave 
term in the expansion in terms of the range of the A^A^ interaction. Thus the present models 
are not perfect but they are good enough for a discussion of the effects of the neutron skin 
and the surface thickness on the proton and neutron scattering data. We compare results 
obtained with three Skyrme interactions SHFl (SKX), SIIF2, and SKM* which will enable 
us to explore the effect of neutron skin (a comparison of SHFl and SHF2) with a fixed 
surface thickness, and surface thickness (a comparison of SHFl and SKM*) with a fixed 
neutron skin. 

We compare also with results obtained with the simple harmonic-oscillator radial wave 
functions for ^''^Pb which were used in [|^. Two sets of oscillator parameters were used. 
For HOI we use hu = 6.70 MeV for both protons and neutrons which is chosen to give the 
rms charge radius of 5.50 fm (a proton rms radius 5.45 fm). HOI has a neutron rms radius 
of 5.84 fm. For H02 the oscillator parameter for neutrons was changed to 7.25 MeV to 
decrease the rms neutron radius to 5.61 fm so that the neutron skin 5* = 0.16 fm is close to 
that obtained with SHFl and SKM*. 

The rms radii from all models of the ground state of ^''^Pb considered are listed in Table |. 
All five give essentially the same radius for the protons but they vary considerably in the 
radius for the neutrons. The difference between the neutron and proton rms radii given in 
the last column emphasizes that spread. Note that the H02 model was chosen to give the 
same rms radii as those of the SHFl model. The neutron radius obtained from the SKM* 
model is similar to those of the SHFl and H02 models. However, as is evident in Fig. |l|, 
each model gives distinctive density distributions. The normalization used is such that their 
volume integrals equate to the proton and neutron numbers of 82 and 126 respectively. In 
Fig. the proton density Pp(r) of both the SHFl and SHF2 models are displayed by the 
solid curve. The density obtained from the SKM* model is given by the double-dot-dashed 
line and exhibits a larger diffuseness compared to the densities of the other Skyrme models. 
Likewise, both the HOI and H02 models have proton densities as given by the dot-dashed 
curve. These quite distinct shapes nevertheless give the same proton rms radius. However, 
they differ in the longitudinal electron scattering form factor as shown in Fig. 0. 
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The five model neutron densities Pni'f') are also shown in Fig. |I[ The SHFl and SHF2 
neutron densities are displayed by the solid and dashed lines respectively. They are similar 
with the SHF2 density having a slightly more diffuse surface region; a property that results 
in the larger neutron rms radius. The density from the SKM* model is given by the double- 
dot-dashed hne and, as for the proton density, exhibits a larger diffuseness compared to 
the densities of the other Skyrme models. The neutron densities of the HOI and H02 
models are displayed by the dotted and dot-dashed lines respectively. As with their proton 
densities, the neutron densities of both of these models are enhanced in the nuclear interior 
over the SHF values. But these HO densities also have increased neutron probability at the 
surface. Recall that the H02 model was set to have the same neutron rms radius as the 
SHFl prescription. 

The SHF (SKX) densities for ^°Ca are displayed in Fig. 0. Also displayed therein are 
the densities from the oscillator model used by Karataglidis and Chadwick [0, for which 
hu = 10.25 MeV. In the surface region, the two models predict essentially the same densities. 
As those densities differ markedly only well within the nuclear volume, we expect influence 
on scattering primarily at high momentum transfer scattering results at energies for which 
absorption through the nucleus is not too large. We anticipate differences in cross sections 
for scattering at energies > 200 MeV and at momentum transfer values > 1 fm~^. 



III. RESULTS 

We have analyzed both proton and (where available) neutron elastic scattering data 
from ^°Ca and ^osp]^ ^g^ gg^ ^^^^^^ 200 MeV. The choices of energies were predicated in 
part on the availability of data and of the momentum transfer values at which those data 
have been measured. In addition our choice was influenced by previous applications of the g- 
folding potentials with those energies being quite successful [0. Furthermore, the effective 



interactions defined for each energy are quite different due to the energy dependence of 
medium effects in the basic g matrices so that the set of NA scattering results we obtain 
provide a consistency check on the various models of structure used. 

The differential cross sections for 40 MeV proton and neutron elastic scattering from 
^"Ca are presented in Fig. |^. Therein, the proton scattering data of Camis et al. 0] and the 
neutron scattering data of de Vito et al. |T9[ are compared to the results of the calculations 
made using SHF (solid line) and HO (dashed line) models. The SHF and HO model results 
equally well describe the data, although they underpredict the proton scattering in the 
regions of the minima. That disagreement may be due to problems in specifying the effective 
interaction at low energies PD|. However, further comment should await the consideration 



of the comparisons of 40 MeV proton scattering from ^''^Pb. Whatever any deficiency at 
this energy there may be though does not affect results we find for, and conclusions we may 
draw from, scattering at 65 and 200 MeV. 

The 65 MeV elastic scattering scattering cross sections for ^°Ca are displayed in Fig. |^. 
Therein, the agreement between the model results and the proton scattering data of Sak- 
aguchi et al. ||2T| is now very good; much better than that found in Fig. ^. Also, there is a 



slight preference for the SHF result at larger angles. In the case of the neutron scattering 



results, both model results agree but underpredict some of the data of Hjort et al. p2 |. 



However, those data are somewhat problematic as concluded from a recent analysis [T^ in 
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which several sets of data were compared with theory at that energy. A new measurement 
of this cross section is required to resolve any such problem. 

It is with 200 MeV scattering from ^°Ca, displayed in Fig. ^ that we observe significant 
differences between the predictions of scattering made using the SHF and HO models. For 
proton scattering that is displayed in Fig. §(a), the SHF model result agrees well with the 
data of Hutcheon et al. [^Sf (circles) and Seifert et al. [0] (squares) and is a significant 
improvement on the result found using the HO model. This variation concurs with our 
expectation that the differences between the inner radial densities of the two models of 
structure would influence the results of scattering calculations at higher energies. Taking 
the results at all three energies, we believe that the SHF model is the better description of 
«Ca. 

We now turn to ^°'^Pb and consider first the integral observables from nucleon scattering 
as a test of sensitivity to the matter distributions of ^°'^Pb. Those integral observables at 
40, 65, and 200 MeV are given in Tables y and ^ and we note that a study of these 
quantities with two of these structure models has been made recently for energies 10 
to 300 MeV. The total reaction cross sections for both proton and neutron scattering from 
^°^Pb are listed in Table 0. Comparing the results of the calculations at both 40 and 65 MeV 
with the available proton data indicates a preference for both SHF, the SKM*, and H02 
models of the ground state density. However, comparing the model results of the neutron 
total reaction cross section with the available evaluated data at those energies p5| gives a 
preference for the Skyrme models only. The predicted proton and neutron total reaction 
cross sections at 200 MeV vary sufficiently that their measurements would be desirable. 

The results of our calculations for the total neutron cross sections are given for the three 
energies in Table IH and are compared with the data of Finlay et al. [2^. For 40 MeV 
the Skyrme models are preferred, although all results overestimate the measured value. At 
65 MeV the Skyrme model results agree well with the data, the SKM* model result doing 
best of all, and do better than both HO predictions. This is not the case at 200 MeV, 
where all results predict the measured value reasonably well. While these total reaction and 
total cross section results together indicate a preference for the Skyrme models, we need 
additional evidence. We consider then the angular distributions of each scattering. 

The differential cross sections for the scattering of 40, 65 and 200 MeV nucleons from 
^°*Pb are presented in Figs. |^, ^ and ^ respectively. At 40 MeV, proton elastic scattering is 
shown in Fig. ^ and evidently there is little if any differentiation between the SHFl, SHF2, 
and SKM* model results. They both compare well with the data of Blumberg et al. P? 



Note, however, that both SHF and the SKM* calculated results agree much better with the 
data than do those found using the HO models of structure. The quality of reproduction 
of the data in this case is in stark contrast to what we found at 40 MeV with ^"Ca. If 
that is to remain a problem with specification of the effective interaction then it seems to 
be a nucleus-dependent effect. In the case of neutron scattering, the results of all model 
calculations agree quite well with the data of de Vito et al. ||19|. The results for ^°^Pb are 
more distinctive than those for ^°Ca at this and other energies, but so are the distinctions 
between the model proton densities for both nuclei. However, these results indicate that 
while the nucleon densities in ^osp^^ ^^^^ better described by the SHF and SKM* models 
nucleon scattering at this energy is largely sensitive to the surface properties only. Only in 
the surface are the proton densities still sufficiently similar for neutron scattering results of 
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all models to be as alike as they are. 

Recall that the integral observables given in Table || make preference to the SHF and 
SKM* models of the density. Similar preference is indicated by the differential cross sections 
for 65 MeV scattering that are displayed in Fig. |^. Considering the neutron scattering first, 
all models give similar results in good agreement with the data of Ibaraki et al. , although 
above 40° the SHF and SKM* model results clearly do better. For proton scattering, the 
SHF and SKM* models are both in good agreement with each other and with the data and 
differ only slightly from those of the HO models. As for 40 MeV scattering, the integral 
observables at 65 MeV concur with these findings that favor the Skyrme models of the 
density. 

The differential cross sections for 200 MeV nucleon scattering from ^°^Pb are shown in 
Fig. ^. In this case only 200 MeV proton scattering data exist for comparison with 
our predictions. Nevertheless, that comparison confirms the findings from analyses of the 
lower energy data, namely that the SHF and SKM* models are the better descriptions of 
the densities of ^°^Pb. In addition, however, some discrimination between the three Skyrme 
model results is evident at this energy. Of the three Skyrme models, the SKM* model 
result agrees best of all with the data while the SHFl does worst. Given that the rms 
radii from these two models are very similar, the marked difference indicates a sensitivity 
to the diffuseness in the density. At this energy, unlike those at lower energies, the neutron 
scattering results show marked differences between the SHF and HO calculations. Again, 
we expect that this is due to the scattering of the higher energy probe being more influenced 
by the bulk nuclear medium properties of the densities. Consequently, a measurement of the 
200 MeV neutron elastic scattering, angular distribution and associated integral observables, 
is certainly desirable. 



IV. CONCLUSIONS 

The distinctions between the predictions of nucleon scattering from ^osp]^ these three 
energies, and found with the set of five model structures used, suffice to select that which 
most likely prescribes the actual matter densities of the nucleus. A similar conclusion is 
reached for ^°Ca for the models considered herein. Specifically we contend that the use of 
(^-folding potential model calculations can differentiate between different model structures 
of the neutron density so as to pin down the neutron rms radius far better than has been 
possible in the past. As well the process gives a good appraisal of the actual density distri- 
bution. For 208p]^ particular, our analyses indicate that the SKM* model gives the best 
representation of the density. Together with analyses of the longitudinal elastic electron 
scattering form factor it suggests a neutron skin thickness for 208p|-, q -^'j f^^. ^ value 
consistent with expectations of the SHFl model nucleus, which is constrained by the FP 
neutron equation of state. The only difference between the two models is that of a larger 
diffuseness for the SKM* model, accounting for both the agreement between the results and 
data for both nucleon and electron scattering. This would suggest a need to extend the SKX 
models to predict a larger diffuseness, for example, by the addition of a (i-wave term in the 
Skyrme Hamiltonian. 

One should also note that while a measurement of the skin thickness as proposed in the 
experiment for the Jefferson Laboratory is important, that quantity is a volume property 
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of the nucleoli distributions. Other information is required to specify a more complete picture 
of the neutron density. From our studies it seems that simultaneous analyses of angular and 
integral observables are relevant. Given that the H02, SHFl, and SKM* models of structure 
we have used predict essentially the same skin thickness for ^°^Pb but give significantly 
different predictions when used to generate NA g folding optical potentials for the nucleon 
scattering, as well as electron scattering form factors, analyses of complementary nucleon 
and electron scattering data permit one to discern such finer details of densities. 
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TABLES 



TABLE L Root-mean-square radii (in fm) for protons (rp) and neutrons 
models are as defined in the text. 



Model 
HOI 
H02 
SHFl 
SHF2 
SKM* 



5.45 
5.45 
5.45 
5.45 
5.45 



5.83 
5.61 
5.61 
5.70 
5.62 



in 208pb. The 



0.38 
0.16 
0.16 
0.25 
0.17 



TABLE II. Total reaction cross sections (in barn) of nucleon scattering from ^''^Pb. The models 



used are as specified in the text. 





40 MeV 




65 MeV 




200 MeV 


Model 


proton 


neutron 


proton 


neutron 


proton 


neutron 


HOI 


2.07 


2.69 


2.11 


2.32 


1.79 


1.77 


H02 


1.95 


2.62 


2.00 


2.27 


1.71 


1.73 


SHFl 


1.89 


2.51 


1.99 


2.19 


1.68 


1.69 


SHF2 


1.95 


2.55 


2.03 


2.22 


1.72 


1.71 


SKM* 


1.92 


2.54 


2.01 


2.21 


1.69 


1.70 


Expt 


2.01 ±0.04 m 


2.50 |25| 


2.02 ± 0.06 M] 


2.20 |25l 







TABLE III. Total cross sections (in barn) of neutron scattering from ^'^^Pb. The models used 
are as specified in the text. 

Model 40 MeV 65 MeV 200 MeV 

HOI 5T0 486 3M 

H02 4.94 4.72 2.97 

SHFl 4.63 4.61 2.94 

SHF2 4.71 4.67 2.94 

SKM* 4.69 4.63 2.96 

Expt |2§ 4.392 ± 0.001 4.634 ± 0.001 2.990 ± 0.003 
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FIGURES 



(a) - 




(b) - 




r (fm) 

FIG. 1. Nucleoli densities in ^''^Pb. The solid and dot-dashed curves in the proton densities pp 
portray both SHF and both HO models respectively. The neutron densities p„ given by the solid, 
dashed, dotted, and dot-dashed lines portray respectively the SHFl, SHF2, HOI, and H02 models. 
The double-dot-dashed line in each case denotes the density obtain from the SKM* model. 
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q (fm ^) 

FIG. 2. Longitudinal elastic electron scattering form factor for ^'^^Pb. The data |31| are com- 
pared to the results of the calculations made using the SHFl, SHF2, and SKM* models portrayed 
by the solid, dashed, and double-dot-dashed lines, respectively. The oscillator result is portrayed 
by the dot-dashed line. 
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FIG. 4. Differential cross sections for 40 MeV nucleon elastic scattering from ^'^Ca. The proton 



scattering data of Camis et al. are compared in (a) with the results of the calculations made 
as defined in the text. In (b), the neutron scattering data of de Vito et al. [|l9j are compared with 
the results of equivalent calculations. 
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FIG. 5. As for Fig. ^ but for 65 MeV scattering. The proton scattering data in (a) are those 
of Sakaguchi et al. ||2l|], while the neutron scattering data in (b) are those of Hjort et al. p^ . 
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FIG. 7. Differential cross sections for 40 MeV nucleon elastic scattering from ^'^^Pb. The proton 



scattering data of Blumberg et al. |27| are compared in (a) to the results from the SHFl and SHF2 
models (solid and dashed lines respectively), the SKM* model (double-dot-dashed line), and to 
the results of the HOI and H02 models (dot-dashed and dotted lines respectively). The neutron 
scattering data of de Vito et al. [19| are compared to the results of those models as defined in (a). 
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FIG. 8. As for Fig. |^ but for 65 MeV scattering. The proton scattering data in (a) are those 
of Sakaguchi et al. ||2l|] , while the neutron scattering data in (b) are those of Ibaraki et al. . 
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